Heavy Quark Fluorescence 
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Heavy hadrons containing heavy quarks (for example, T-mesons) feature a scale separation be- 
tween the heavy quark mass (about 4.5 GeV for the &-quark) and the QCD scale (about 0.3 GeV) 
that controls effective masses of lighter constituents. Therefore, as in ordinary molecules, the deex- 
citation of the lighter, faster degrees of freedom leaves the velocity distribution of the heavy quarks 
unchanged, populating the available decay channels in qualitatively predictable ways. Automatically 
an application of the Franck-Condon principle of molecular physics explains several puzzling results 
of T(5S') decays as measured by the Belle collaboration, such as the high rate of B*Bg versus B*Bs 
production, the strength of three-body B*Btt decays, or the dip in B momentum shown in these 
decays. We argue that the data is showing the first Sturm-Liouville zero of the T{5S) quantum 
mechanical squared wavefunction, and providing evidence for a largely bb composition of this meson. 



The Franck-Condon principle [l| is essential in under- 
standing fluorescence in the decay of excited molecules. 
It states that, during an electronic transition, the elec- 
tron orbital relaxes to the ground state in a time too short 
for the nuclei to react. Hence the wavefunction associated 
to the motion of the nuclei remains the same after the 
transition. Since the ground-state wavefunctions of this 
motion in the two electronic adiabatic potentials (before 
and after the transition) are unequal, the nuclei remain in 
a superposition of excited-states, resulting in secondary 
transitions (efficiently degrading excitation energy). 

While the Born-Oppcnhcimer approximation from 
molecular physics had already been studied for decades 
in the context of particle physics, we have only recently 
understood Q the relevance of the Franck-Condon prin- 
ciple here. The key is its usefulness to analyze open flavor 
decays of heavy hadrons. Such decays are characterized 
by the heavy quarks from the initial hadron streaming 
away in separate hadrons. Wc take as example the newly 
analyzed T(10860) ^ BB, BBtt . . . decays. 

The momentum of the heavy quark inside the heavy 
hadron, fcf,, follows a quantum- mechanical probability 
distribution that depends on the wavefunction of the 
state. This wavefunction accepts a decomposition in 
an appropriately chosen quark-gluon basis, for example 
|T) = ao\bb) + ai\bbg) + a2\bbqq) + a^lggg) + ... 
One of the salient problems of hadron physics nowadays 
is to determine the functions a^, whose size character- 
izes whether a state is "mostly a quarkonium" or other- 
wise. Characteristic of pure (excited) quarkonia is the 
presence of Sturm-Liouville zeroes in the momentum dis- 
tribution. This follows from the Hermiticity of the QCD 
Hamiltonian that implies orthogonality between dif- 
ferent wavefunctions. Fig. [T] is a realization of T{5S) 
wavefunction (candidate bb assignment for T(10860)) in 
a model approximation to Coulomb gauge QCD One 
clearly sees the 4 zeroes of the 55 excitation, beautifully 
exploited in other areas such as nuclear physics @. 

The advantage of open-flavor decays (experimentally 
more difficult to reconstruct than transitions to ground- 
state T) is that the velocity Vb of the initial-state quark 
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FIG. 1: Top: Coulomb-gauge model wavefunction of T(55'). 
Bottom: Data points are 2-body decay widths of the T( 10860) 
(Belle collaboration), plotted as a function of the center of 
mass momentum for the decay products. From right to left, 
BB, BB* + B*B, B'B*, BsB^, BsB* + B^Bs, We 
have divided the experimental data by known spin, isospin 
and phase-space factors so it can be directly compared to 
the square of the top panel wavefunction (solid line) or to the 
wavefunction folded with a complete ^Pq decay model (dashed 
line), that gives us an idea of the errors of order Aqco /rrii, in- 
volved in a naive application of the Franck-Condon principle. 



approximately equals vb, the velocity of the final state 
meson. This is because the final i3-meson momentum 
k in a transition driven by gluons and light quarks sat- 
isfies k = kf, + ©(Aqcd), which is the rationale behind 
effective-theory analysis such as Heavy Quark Effective 
Theory (HQET) di]. 



2 



The 2-body data T(55') ^ consists of the 

partial 2-body decay widths. To obtain them from the 
decay ampUtude one integrates over phase space 



dT = 



(27r)4(5W(P-/fc-/fc') 
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where P, k and k' are the 4-momenta of the T, B and 
B, respectively. Integration yields the well-known 



r(fc) 



}f_ 

Stt Mr EkEk' 



\M\' 



(2) 



where the modulus of the B-mcson momentum k (in 
the CMS) is fixed by kinematics to Kallen's function 
k = j^X^/^{M^,A''P,M'^), with M,M' either of 
Mb, Mb* ,Mb,, Mb'- Thus each 2-body decay probes 
the kb distribution in the initial state only at one point 
(blurred by O(Aqcd)). But since there are six al- 
lowed decay channels, we have access to six different 
momenta [9|, allowing the suggestive mapping of the 
T(10860) fcb distribution from experimental data in Fig.[T] 
(ignore for now the calculated curves). Automatically 
one understands the puzzling feature that the decay 
width to B*B* (19 MeV) is a factor 12 times the width 
to BsB* + B^Bs (1.6 MeV), instead of the (7/4) ratio 
expected from spin counting alone [loj . We can now un- 
derstand that the second decay just happens to fall near 
the first Sturm-Liouville zero of the 5S'-wavefunction in 
the initial state, whereas the first is well into the wave- 
function maximum. Thus the anomalously large strange- 
vector-vector decay finds a simple explanation in terms 
of the structure of the parent hadron. 

If one wishes better resolution of the kf, distribution 
inside T(10860), resort needs to be made to 3-body de- 



cays, from which B^*'^B'^*\, first analyzed in [11|, is the 
more promising given that T(55) is not so high above 
open-flavor threshold. In these decays the pion carries 
off a momentum p that balances momentum conserva- 
tion. Then if k' (the B^*^ momentum) is constrained by 
energy conservation, k (the B*^*^ momentum) is a free 
variable within its phase-space bounds. The partial de- 
cay width as a function of this k then allows access to k^ 
in the initial state. Performing immediate angular inte- 
grals, and leaving k ^ kf, free, it is easy to show that 
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dk 
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where fc^^^ is the maximum value of the i3-meson mo- 
mentum given by the kinematics and x' is the polar angle 
cosine of k', fixed by energy conservation to 
2kk'x' = {Mr - Ek - Ek'f - (M^ + k^ + k'^). If 
the T(5S') has a significant bb component, the decay 
can be interpreted at the quark level as bb ^ (6 + 
light) (6 + light) (light). Then the Franck-Condon prin- 
ciple (up to corrections Aqcd that we can estimate 




FIG. 2: Dashed line: computed momentum distribution 
of the B^*'-meson in three-body T(5S') B^*^B'-*^tt non- 
strange decays (summed over all B/B* spins). The vari- 
able Mbc is related with the B meson momentum as Mtc = 
\/ (A/t/2)2 — fc^.The momentum of the 6 quark in the parent 
T(55) is given hy kb — k + O(Aqcd). Thus, except at very 
low k, this distribution probes the confined quark degrees of 
freedom in the initial state. The dip in k at around 600 MeV 
(or Mbc = 5.4 GeV in the Figure) is caused by the first Sturm- 
Liouville zero in the T(55') kb distribution (solid) folded with 
3-body phase space and is correlated with the small BsBl 
decay width (see Fig. [T]). This zero is indeed visible in the 
recent data (Fig. 2b from [l^'). 



within a model) states that the velocity distribution of 
the final-state B meson is the same as the velocity dis- 
tribution of the 6 quark in the initial T(5S') state, thus 
providing us with a window to its structure. In Fig. [2] 
we show the calculated momentum distribution of the 
B meson in these decays (this is directly measurable by 
Belle; to obtain a velocity distribution by eye inspec- 
tion, note vb — Ps/'niB is quite a good approximation 
at these energies). A dip is clearly visible at around 
Mbc = 5.4 GeV (or k = 600 MeV). This is a consequence 
of the first Sturm-Liouville zero in the T{bS) wavefunc- 
tion and is correlated with the small B^B* decay width, 
since the momentum for that 2-body decay is also around 
600 MeV. Remarkably, structural insight from the parent 
hadron, in terms of confined degrees of freedom, helps to 
interpret puzzling features of decay data. Only the lowest 
kb < AqcD is not accessible (application of the Franck- 
Condon principle just shows that the heavy quarks re- 
main at rest up to O (Aqcd/ 'm-b)). 

In terms of the Franck-Condon principle, the reported 
large fraction of 3-body decays of T (10860) 0], that 
has not been quantitatively predicted by theoretical ap- 
proaches [l^l, is easy to understand. Indeed, in these 
decays there are always final-state configurations avail- 
able whatever the kb in the initial state might be, since 
the pion momentum is an additional free variable to en- 
sure conservation laws. (In 2-body decays on the other 
hand, only one final k is open, so that only small parts 
of the kb wavcfunction overlap with the decay operator). 

We now present some calculational detail for the model 
theory curves provided to illustrate the principle. 
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The structure of the initial T(5S') and final states B, B*, 
IT, etc. is taken from the Cornell model of Coulomb-gauge 
QCD based on the Hamiltonian for the quark fields 



rfx*^(2;)a-A(x)*(x) 
rfx*;f(x)(-ia- V + /3mJ*g(x) 
dxdyp''(x)V^(|x-y|K(y) . 



(4) 



The a ■ A term in the first line is reduced to a trans- 
verse potential between quark sources (of Yukawa type) 
as explained in Q and the Cornell (Coulomb+linear po- 
tential) is used for the Coulomb static kernel V, a model 
philosophy consistent with lattice [isj results for very 
heavy quarks. This relativistic many-body Hamiltonian 
is treated in BCS+TDA approximations to yield the 
dressed bq or bb canonical wavefunctions for pseudoscalar 



Xl^{\ci\)=Yoo{ci) 



and vector mesons (with spin index a 
ll) = >oo(q) 
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(We have also computed the D-wave component and 
found it to be totally negligible in bottomonium). We 
further neglect the small model difference between the 
radial B and B* wavefunctions (which is again supported 
by HQET for large mt). 

The values of the model parameters are now given, gluon 
mass nig = 0.89 GeV (string tension, a is dependent 
as per the cited work of Q), strong Colomb coupling 
gs = 1.41 and heavy-quark mass mt, = 3.73 GeV. These 
parameters yield good agreement with experimental data 
on the T bb spectrum as sketched in Figure [31 
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FIG. 3: T spectrum in Cornell-type model of Coulomb-gauge 
QCD, in agreement with the experimental data [13 ]. 

As for the decay mechanism, a complete QCD- 
rescattering calculation is beyond our ability, so we sup- 



plement the model with the time-tested ^Pq prescrip- 
tion [l6j . In this framework, light quark creation pro- 
ceeds via a chiral-symmetry breaking interaction Hj pro- 
portional to (T • p. Dyson-Schwinger studies [l8| have 
shown that the quark-gluon vertex features strong scalar 
contributions only up to momenta of order Aqcd , where 
the light quark running mass function M{k) falls steeply 
like a power-law towards the current-quark mass mq. 
Thus we take the ^Pq vertex proportional to the BCS 
mass-gap M{k)/E{k) generated by Eq. that signals 
spontaneous chiral symmetry breaking, with a propor- 
tionality constant c serving as decay-model parameter. 

The Feynman amplitude ([mass dimension] =1 for the 
2-body decay) parametrizes the S-matrix (5=1-1- zT), 



{kk'\T\P) = (27r)4(5W(P - fc - k')M 



(7) 



The Born scries iT = —i J dtHj + • 
the ^Pq decay-model Hamiltonian 



connects it then to 



dt{kk'\Hi\P) = (27r)4<5('*)(P -k- k')M (8) 
(Pq M{q) 



(2^)3 E{q) 
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The time independence of Hj leaves, as customary, the 
decay rate limr^^ J^j,{kk'\Hi\P)dt = S{0){kk'\Hi\P) . 
The decay amplitude is then depicted in Fig. [H (left). Its 
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FIG. 4: Model T 2- and 3-body open-flavor decays through 
^Po vertices. In the limit Aqcd/A/(, — they respect the 
Franck-Condon principle that QCD should also obey, and the 
momenta of the heavy quark lines flowing from the left is close 
to the external momenta of the decay mesons on the right. 



computation requires the matrix element in Ec^. ^ 

M{\k-q\ 



(fcfc'|Pj|P) = c<5(3)(k + k') J rf^q. 



P(|k-q|: 
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The Franck-Condon principle is patent in this ex- 
pression; if the decay vertex M(|k — q|) is small for 
|k — q| » Aqcd — 250 MeV, the heavy-quark momen- 
tum k — 2q, which is twice the argument of the P-meson 
wavefunction in Eq. (jlOp . becomes the external me- 
son momentum k up to O(Aqcd)- The heavy-quark and 
meson velocities are then approximately the same as ex- 
pected. Thus we have a simple model realization of the 
general principle. 
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We average over initial and sum over final spin states 
for unpolarized decays. Particularizing to BB, we have 



(11) 



In 2-body channels the spin factorizes from the integral 
for all three final state spin combinations, leaving the 
known ratios 1:4:7 for BB : BB* + B*B : B* B* . 

The 3-body decays depicted on the right of Fig. U] are 
an immediate extension of these formulae, except that 
the pion wavefunction requires an additional word of cau- 
tion. It is computed in the rest frame defined by the 
Hamiltonian in Eq. Q but the velocity of the outgoing 
pion is relativistic. While rigorous studies of the boost 

this 
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operator in Coulomb QCD are underway 
work we have limited ourselves to a kinematic Lorentz 
contraction of the computed wavefunction. 

Looking forward to further refinements of this work, we 
note that our estimate of the divergence from the Franck- 
Condon principle is model-based, and perhaps one could 
directly use effective theory language to parametrize this 
irreducible error. But it is not currently clear what effec- 
tive theory one should pursue. 

Ground state hh is deeply bound respect to the open- 
flavor decay threshold. Typical momentum transfers are 
of order QsMf, and do not vanish asymptotically for large 
quark masses at fixed , the appropriate effective theory 
being NRQCD. However, for resonances above threshold 
such as T(10860) some extension of HQET might ap- 
ply, since the open flavor B mesons perceive momentum 
transfers from the vertex to the heavy quark of order 
the string tension scale, not proportional to m^. This 
deserves further study. In both cases though, the lead- 
ing order phenomenon is the equality of heavy quark and 
heavy meson velocity that we have exploited. 

On the phenomenological side, a large number of 
new resonances in the high-charmonium spectrum, 3.8 — 
4.6 GeV, has been unveiled in the last years. These 
new states, while pending classification named X(3940), 
r(4260), r(4320), ... or generically X - Y - Z reso- 
nances, are seen either recoiling against & J/ip (such as 
the 3940) or decaying to J/i> + mr with n = 1, 2, 3 . . . 

Unstable resonances pose very challenging problems to 
lattice gauge theory. No known effective theory can de- 
scribe either of all these resonances from first principles 
without additional model assumptions. The amount of 
"hard" theoretical understanding is then very limited, 
and the discussion is still at the level of simple models. 

Our observation is that some structural information 
(such as momentum distributions of heavy quarks inside 
the new states) can be obtained from open-flavor decays. 
Several examples [1] have been presented of the possible 
use of the fluorescence analogy to analyze the decays of 
the new charmonium-like states XYZ. A (presumedly) 



excited cc state would feature Sturm-Liouville zeroes in 
the decay ip* — > DDir, by tracking the momentum of the 
D. An exotic meson (hybrid, tetraquark, etc) of equal 
mass would not present such Sturm-Liouville nodes in 
its relative cc momentum distribution. 
Another application consists of constructing an off- 
planes correlator in four-body decays to identify the like- 
liness of a tetraquark component in the excited meson. 
Yet a third idea is to use DD momentum distributions 
in certain B-meson decays to learn about the charm-sea 
in hadrons. We look forward to turning some of these 
possibilities into usable experimental tests. 
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